Ion Landau Damping on Drift Tearing Modes 



00 



Of 



J W Connor 1,2 ' 3 , R J Hastie 1 ' 3 and A Zocco 1 ' 3 
1 EURATOM/CCFE Fusion Association, 
Culham Science Centre, Abingdon, Oxon, 0X14 3DB, UK 
2 Imperial College of Science Technology and Medicine, London, SW7 2BZ, UK 
3 Rudolf Peierls Centre for Theoretical Physics, 
IKeble Road, Oxford, 0X1 3NP, UK 

The equations governing the ion Landau damping (ILD) layers for a drift tearing 
mode are derived and solved to provide a matching to ideal MHD solutions at large 
x and to the drift tearing solution emerging from the ion kinetic region, kpi ~ 1, at 
small x, the distance from the rational surface. The ILD layers lie on either side of the 
mode rational surface at locations defined by k y xVri/L s = £J* e (l + 0.73r/ e ) and have 
been ignored in many previous analyses of linear drift tearing stability The effect 



of the ILD layer on the drift tearing mode is to introduce an additional stabilizing 



contribution, requiring even larger values of the stability index, A' for instability, 
than predicted in refs.[l|, 2]. The magnitude and scaling of the new stabilizing effect 
in slab geometry is discussed. 

I Introduction 

Tearing modes are driven by current gradients and pressure gradients in the plasma, 
whose destabilising effect is characterized by a quantity A'. The associated tearing of the 
magnetic field occurs at the resonant surface, ro, where m = ng(ro), m and n being the 
poloidal and toroidal mode numbers of the helical instability. In the vicinity of r = ro 
complex kinetic effects can occur. To account for these it is convenient to treat a narrow 
'inner' region in which these effects are included but the equations are simplified by virtue 
of the localised nature of this layer. Stability is then determined by matching solutions 
of these inner equations to outer ones, valid across the regions away from the layer (and 
satisfying appropriate boundary conditions at the plasma boundary and magnetic axis) 
where a simple ideal MHD description of the plasma suffices. The solution of the inner 



2 



equations provides a quantity A (a;), where u is the mode frequency, while the solution of 
the outer equations yields the, aforementioned, quantity A': the matching condition leads 
to a dispersion relation. 

For a simple resistive MHD model, (or even a cold ion model containing electron 
diamagnetic effects) of the inner region, the corresponding solutions match satisfactorily 
to the outer ones, with the perturbed parallel electric field, En vanishing as one reaches 
the outer ideal MHD region. However, for a kinetic treatment of ion physics, when finite 
ion Larmor radius, FLR, and parallel thermal motion are included, this matching becomes 
more complex. Even though one can match the functional dependence of the longitudinal 
component of the perturbed vector potential, An from the inner region to its form in the 
ideal MHD region, where An ~ 1 + (r — r )A', the inner region expression for E\\ fails to 
vanish as it should. In order to accomplish this (i.e. decay of En ) within the matching 
range, one must consider an intermediate region where ion sound and ion Landau damping 
become significant at the characteristic frequency of a drift tearing mode. As a consequence 
the overall matching condition is modified with an impact on the stability of the drift 
tearing mode. The investigation of this effect in a cylindrical plasma model is the purpose 
of this paper. 



II Plasma Model 

n 

We start from the fourth order coupled set of equations[l], in slab geometry, comprising 
the vorticity equation: 

feu J {{ = [(i - r (6))(w - oo») - u>tfnb(r - ro] , (i) 



and Ampere's equation, 



d 2 A 



-~J\\, (2) 



dx 2 

where and A\\ are the perturbed electrostatic potential and longitudinal component of the 
perturbed vector potential, fey = k y x/ L s is the longitudinal wavenumber, with x the distance 
from the mode rational surface and L s shear length: in toroidal geometry L s = Rq/s, with 
R the major radius, s = rq'/q and q the safety factor. In eqs.1,2, Jn = Jiu + J\\ e is the 
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perturbed current density, e, the proton charge, n, the equilibrium density, u the mode fre- 
quency, = k y Tj/(ejBL n ) the diamagnetic frequency of each species, 1/L n = dLog (n) / dx , 
rjj = dLog(Tj) / dLog(n); T n (b) = I n (b)e~ b , b = k^Ti / (mju;^) , with I n the modified Bessel 
functions; k\ = k 2 + k 2 , denotes the component of the wave number perpendicular to the 
magnetic field, with k y the component of k± within the magnetic surface, and k the compo- 
nent perpendicular to the surface. Finally, u ci denotes the ion cyclotron frequency. 
Expressions for the perturbed current densities, J^ e and Jn ti are derived from the respective 
kinetic equations for electrons and ions (in the kpi <C 1 limit): 

e 2 nE\\ , 



■h 



— I- 



klT P 



(3) 



J\\A = 



e 2 nE\\ 



- UH r ]i {C 2 {l + C i Z{C i ))--C i Z{Q) 



(4) 



where Q = ooLs/kyxVxi, Vt 4 = \JlTijmi, Z is the plasma dispersion function. Making 
use of the long wavelength approximation for the ions in eq.l, with b ~ p 2 d 2 /dx 2 <C 1, 
eliminating J\\ i e and introducing a dimensionless variable X = we obtain the following 
fourth order system describing the potentials 0, An, which we now denote by A.: 



X 
X 



d 2 A 
dX 2 
d 2 A 
dX 2 



L 



u*i d 2 (f) 
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u 1 dX 2 ' 

2 

s W 



(5) 
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UJ 
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1+ l 



w*i rji f Z X Z 

- T -x- 2 \ l + x- — 



(6) 



where = (f)/V T i, Z = Z(l/X), r = T e /Ti and f3j are the ratios of plasma to magnetic 
pressure for each species. 



Ill Solving the fourth order system through the Ion Landau region. 

To solve the fourth order system of eqs.5 and 6 we first exploit the small parameter Pi ^ 1 
in eq.5 to obtain; 
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which matches correctly to the ideal MHD outer solution as X — > oo. Then, inserting this 
result in the RHS of eq.6, we use eq.5 to generate the following second order inhomogeneous 
ODE for 0: 

Q{X) 



dx 2 



1 

X 



A' ojL, 



2 kyVn 



where, 

Q(X) 



OJ 



2L 2 n uT + l + rji 

OJ 

0J*p. 



Cj-l+{ujT + l){l + —Z)+r H 



1 1 7 Z 
X 2 ~ + X 1 ~2X 



(8) 



(9) 



and we have dropped the caret on for simplicity. 

The objective is now to solve eq.8 throughout the ion Landau region where X ~ 1 
in order to determine the asymptotic form of <fi at small values of X, corresponding to 
Pi <C x <C pi(L s /L n ). In this region <p(X) will match to the drift tearing solution emerging 
from the ion-kinetic region, where x ~ pj. This also takes the form <ft ~ (a + 1/x) and in 
most earlier analyses of drift tearing modes, e.g. [2| , a = A'/2 has been assumed; however 
this requires justification. For this purpose we assume that the relevant mode frequency is 
known and is characteristic of the drift tearing mode in the semi-collisional electron limit 
i.e. oj = u me (l + 0.73r] e ). 

Noting that Q(X) ~ 0{r 2 L 2 s /L 2 n ) > 1 throughout the region, we solve eq.8 in the WKB 
approximation so that the full solution is constructed from the two WKB solutions of the 
homogeneous eqn.: 

1 



MX) 



Q l/i 



■ /(f a/ Q(y)dy p - / X yf Q(y)dy 



(10) 



The relevant solution, satisfying the matching condition to ideal MHD for X ^> 1, takes the 
form: 



0(X) 

r(v) = Q(v) 



l -{n (Y) / r(!j)n_ {! j),l!j + 0-{X) 



A' 



' x 

ojL s 1 

2kyV Ti + X 



(11) 

(12) 



In eq.ll the factor 1/2 comes from the Wronskian, W((j)+,<fi-) = —2 and c and Xq are 
arbitrary constants. In practice it will be convenient to choose a value of X Q 1 when 
evaluating the asymptotic form of (f)(X) in the X <^ 1 region, in the next section. 
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IV Asymptotic form of <fi in the small X limit. 

In order to match to the drift tearing mode solution in the region x ^> Pi we require the 
functional form of <f){X) in the region \j\fQ <I<1. In this limit, the ion Z-functions 
can be approximated by the "ion-sound" expansion with the result: 

+ l + d3) 

= Qo [1 - AX 2 ] , (14) 
1 f utL s \ 2 u-l lc2;r + l+?7j 

Q ° ~ 2{—J ur+l + r,^ 2 u-1 (15) 

Now, making use of these results within eq.ll, we find that 0_(X) ~ e -VQo x j s exponentially 
small so that the second term in eq.ll (containing the constant c and the integration end 
point Xq ) can be neglected in the evaluation of 0. The integral in the first term of eq.ll is 
evaluated by successive partial integrations, leading to the result: 



4>{X) 



A' ujL, A 



4, (16) 



2 kyVn y/Qo 

to first order in the small parameter, 1/y/Qo ~ L n /L s 

We translate this result back into an "effective" value of A' to be used for the matching 
of a drift tearing solution, to an outer (in x) solution. When this matching is performed 
within the ion Landau region, one finds this effective A' is given by: 

L n \ 2 f cut + 1 + rji \ 3/2 2 



*4,* = "»-{t.) I^t^J (17) 

where, as noted earlier, Cj — 1 + 0.73?] e for a drift tearing mode. 



V Discussion and Conclusions 

The result obtained in eq.17 implies that the ion sound/ion Landau region introduces a 
stabilising effect for drift tearing modes. This is in addition to the diamagnetic stabilisation 

n n 

found by Drake et al.[3] and the finite ion orbit stabilisation found by Cowley et al.[l| and 



a. 



in [2]. The resulting stability threshold for drift tearing can now be expressed in the form: 
A'pi > A FLR + A DIA + A ILD , (18) 



(u-1) 2 (lu + l)(w + 1-^/2) 



A FLR = ^ log{Y)-- , (19) 
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A DIA = n$(u-l)\I\, (20) 
2 /w + l + r/A 3/2 2 



A '- = W {^r 1 ) w (21) 

where F ~ {Pi/$) ^ 1 with 5 the width of the electron current layer, 7(d), r/j) is an integral 
defined in eq.24 of [2] (see also Fig.4 of Q), /3 = (/3 e /2)(L^/L^) and for simplicity we have 
set r = 1. The expressions in eqs.19 and 20 apply when the parameter, /3 < 1. At high 
values of /9, complete screening of the resonant surface occurs resulting in elimination of the 
drift tearing mode. We conclude that the additional stabilization represented by Ajld can 
only be of significance at very small values of /3 e or at small values of T] e . 
Similar "ion sound" modifications of the effective A' for drift tearing modes have previously 
been derived by Bussac et al.Q] and by Fitzpatrick[5] using fluid treatments. In ref.^|, semi- 
collisional drift tearing modes were considered in the small Larmor radius limit, 5. Slab 
geometry was assumed and temperature gradients were (initially) neglected. The resulting 
expression for the stability threshold takes the form: 



A' Pi ~ & 



1/2 



(22) 



and it was argued that the stabilizing effect becomes much weaker in the presence of an 
electron temperature gradient, i] e ^ 0. This same trend is apparent in eq.21 above, where 
the limit rj e = appears to be singular. This limit would require a calculation of 5u, a small 
correction to the drift tearing frequency, Uq = 1, as was carried out in ref.J^I in the pi ^ 5 
limit. We do not pursue this special case, where 7] e = 0, but note that this explains the 
differing results in eqs.21 and 22. In j5| toroidal effects were included and analytic solutions 
generated in a transition region, between the electron current channel and the MHD region, 
by introducing, as a simple model, a Heaviside function switch-on of the perturbed ion 
current. Using this device a result of the form: 

A' eff = aA' + 5, (23) 

was obtained, where the multiplier, a is positive, greater than unity, but reduces to unity 
in the slab limit. The offset term, 5, found by Fitzpatrick is positive (destabilising), but 
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vanishes in the slab limit. 

The foregoing analysis differs from both these treatments in being a kinetic treatment 
including finite ion orbit effects and in which ion Landau damping can occur. However, the 
modification to A' generated by the ion transition region is not complex; i.e. the imaginary 
part of the ion Z function never features in the above analysis. The apparent cause of this 
can be identified by investigating the variation of the longitudinal electric field, E\\ through 
the ion transition region. It appears that E\\, which must tend to zero as the outer ideal 
MHD region is approached, actually becomes small as the ion sound region is traversed by 
the drift tearing eigenfunction, i.e. before the ion Landau resonances are encountered. This 
indicates that a fluid treatment of the ion transition region should suffice. 
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